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Abstract: Snow cover is an important water source for vegetation growth in arid and semi-arid areas, and 
grassland phenology provides valuable information on the response of terrestrial ecosystems to climate 
change. The Mongolian Plateau features both abundant snow cover resources and typical grassland 
ecosystems. In recent years, with the intensification of global climate change, the snow cover on the 
Mongolian Plateau has changed correspondingly, with resulting effects on vegetation growth. In this study, 
using MOD10A1 snow cover data and MOD13A1 Normalized Difference Vegetation Index (NDVI) data 
combined with remote sensing (RS) and geographic information system (GIS) techniques, we analyzed the 
spatiotemporal changes in snow cover and grassland phenology on the Mongolian Plateau from 2001 to 
2018. The correlation analysis and grey relation analysis were used to determine the influence of snow 
cover parameters (snow cover fraction (SCF), snow cover duration (SCD), snow cover onset date (SCOD), 
and snow cover end date (SCED)) on different types of grassland vegetation. The results showed wide 
snow cover areas, an early start time, a late end time, and a long duration of snow cover over the northern 
Mongolian Plateau. Additionally, a late start, an early end, and a short duration were observed for 
grassland phenology, but the southern area showed the opposite trend. The SCF decreased at an annual 
rate of 0.33%. The SCD was shortened at an annual rate of 0.57 d. The SCOD and SCED in more than 
half of the study area advanced at annual rates of 5.33 and 5.74 DOY (day of year), respectively. For 
grassland phenology, the start of the growing season (SOS) advanced at an annual rate of 0.03 DOY, the 
end of the growing season (EOS) was delayed at an annual rate of 0.14 DOY, and the length of the 
growing season (LOS) was prolonged at an annual rate of 0.17 d. The SCF, SCD, and SCED in the snow 
season were significantly positively correlated with the SOS and negatively correlated with the EOS and 
LOS. The SCOD was significantly negatively correlated with the SOS and positively correlated with the 
EOS and LOS. The SCD and SCF can directly affect the SOS of grassland vegetation, while the EOS and 
LOS were obviously influenced by the SCOD and SCED. This study provides a scientific basis for 
exploring the response trends of alpine vegetation to global climate change. 
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1 Introduction 


Snow cover is a unique natural attribute of the cold season at middle and high latitudes; it is also 
an important part of the cryosphere on Earth and can affect the surface water supply, global 
energy cycle, and atmosphere-freezing cycle (Liston, 1999; Han et al., 2014). Affected by climate 
change, the global snow cover has changed significantly, and the snow cover area in the Northern 
Hemisphere has shown a decreasing trend (IPCC, 2008; Brown and Mote, 2009; Brown and 
Robinson, 2011), which not only affects soil water, heat, and the radiation energy balance directly 
but also indirectly influences the ecosystem balance. Vegetation is a critical part of the ecosystem, 
and its phenological changes are regarded as a comprehensive "indicator" of ecological 
environmental change (Hou et al., 2014a; Kong et al., 2017), which has become the core issue of 
global change and carbon cycle research (Piao et al., 2006; Shen et al., 2014). Vegetation 
phenological changes are generally believed to be closely related to climate change. Moreover, in 
recent years, the experimental and observational results of many scholars have shown that 
vegetation phenological changes are closely correlated with seasonal snow cover (Thompson et 
al., 2015; Jin et al., 2017), and winter snow cover changes also affect plant growth and 
productivity (Peng et al., 2010; Wipf and Rixen, 2010; Matsumura and Yamazaki, 2012; Trujillo 
et al., 2012; Mark et al., 2015; Thackeray et al., 2016; Wang et al., 2018a). 

Snow cover is an important water resource in arid and semi-arid areas. Important snow cover 
parameters, such as snow cover season (SCS), snow cover fraction (SCF), snow cover duration 
(SCD), snow cover onset date (SCOD), and snow cover end date (SCED), can reflect the climatic, 
hydrological, and ecological conditions of a certain area directly or indirectly (Thackeray et al., 
2016; Wang et al., 2016). Based on the 8 d snow cover products derived from MODIS (Moderate 
Resolution Imaging Spectroradiometer) Terra/Aqua data, Liu et al. (2017) found that the 
inter-annual variations of snow cover presented a slight and non-significant increase trend from 
2000 to 2015 for the whole Central Asia. Chen et al. (2016) pointed out that in China, the cooling 
effect of snow cover was enhanced from 1983 to 2010, which indicates that snow cover does 
impact climate. Simultaneously, the thaw time of soil and its subsequent soil moisture content 
were found to be the main influencing factors of vegetation phenology in temperate and cold 
zones (Shen et al., 2011). After a long and cold winter, the warming temperature and meltwater 
promote vegetation growth, especially in alpine region (Dorji et al., 2013). Peng et al. (2010) 
noted that in China, winter snow cover plays a crucial role in regulating desert vegetation growth, 
and there is also a significant positive correlation between winter snow depth and grassland 
Normalized Difference Vegetation Index (NDVI). Wang et al. (2018b) found that snow cover 
phenology has a strong controlling effect on the growth dynamics of alpine vegetation: the longer 
the SCD of alpine grassland lasts, the earlier the start of growing season (SOS) appears, and the 
longer the length of growing season (LOS) lasts. In addition, snow cover has a thermo protective 
effect, ensuring that vegetation can survive in harsh winter environments. 

The Mongolian Plateau with abundant snow cover resources and typical grassland ecosystems 
is in the middle of Asia and belongs to the arid and semi-arid climate zone. The Inner Mongolian 
Plateau is the most important green ecological barrier and one of the snow-rich regions on the 
Mongolian Plateau (Sa et al., 2012; Bao et al., 2014; Dong et al., 2019). Snow cover varies 
regionally and has different effects on different vegetation types (Shen et al., 2016; Qiao and 
Wang, 2019). There are few reports on snow cover on the Mongolian Plateau and its relationship 
with grassland phenology. 

In this paper, we extracted snow cover and vegetation phenological parameters based on 
MOD10A1 snow data and MOD13A1 NDVI data, respectively, to reveal their spatiotemporal 
variation characteristics on the Mongolian Plateau. Correlation analysis and grey relation analysis 
(GRA) were used to analyze the potential influence of the spatiotemporal variation of snow cover 
parameters on grassland phenology and to explore the possible reasons for the relationship 
between snow cover and grassland phenology. The results of this study are of considerable 
relevance for the long-term ecological protection of the Mongolian Plateau and the formulation of 
countermeasures to adapt to global climate change. 
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2 Materials and methods 


2. Study area 


The Mongolian Plateau generally refers to the inland plateau of Northeast Asia, from the Da 
Hinggan Mountains in the east to the Altay Mountains in the west, and from the Sayan Mountains 
and Yablonoi Mountains in the north to the Yinshan Mountains in the south (37°22'-53°20'N, 
87°43'-126°04'E; Fig. 1a). The Mongolian Plateau comprises Mongolia, parts of southern Russia, 
and northern parts of Inner Mongolia Autonomous Region and parts of Xinjiang Uygur 
Autonomous Region of China. The terrain of the Mongolian Plateau is high in the west and low in 
the east, with an average altitude of 1580 m asl. The average annual precipitation is 
approximately 200 mm. The annual average temperature range is large and primarily above 40?C. 
The highest temperature in summer is over 30?C, and the lowest temperature in winter is -45?C. 
The natural vegetation from northeast to southwest is forests, grasslands, and deserts. Grassland, 
as a typical natural vegetation community of the Mongolian Plateau, can be divided into meadow 
steppe, typical steppe, desert steppe, and alpine steppe (Bao et al., 2013). The distribution of 
grassland types on the Mongolian Plateau is shown in Figure 1b. 
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Fig. 1 Overview of the Mongolian Plateau (a) and distribution of grassland types (b) 
2.2 Datasets 


Daily snow cover data at a spatial resolution of 500 m were obtained from the MOD10A1 
Collection 6 (C6) snow product derived from the National Snow and Ice Data Center (NSDIC; 
https://nsidc.org/data/MOD 10A 1/versions/6), with a time range from September of the current year 
to April of the next year during 2000-2018. The MOD10A1 data were batched using the MODIS 
Reprojection Tool (MRT). The geographic coordinates were converted to World Geodetic System 
84 datum, the projection was converted to Albers projection, and the format was converted to 
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Geo-TIFF format, and then, the image mosaic was completed. 

The MODIS NDVI dataset is based on the MODIS NDVI product of MOD13A1 C6 with the 
temporal resolution of 16 d and spatial resolution of 500 m; it was obtained from the Application for 
Extracting and Exploring Analysis Ready Samples (AppEEARS) by Land Processes Distributed 
Active Archive Center (LP DAAC) at the United States Geological Survey (USGS), with a time 
range from 2001 to 2018. Vegetation type classification (meadow steppe, typical steppe, desert 
steppe, and alpine steppe) data of the Mongolian Plateau were obtained from Bao et al. (2017). 

The Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) Global 
Digital Elevation Model 2 (GDEM2) derived from the NASA's LP DAAC 
(https://Ipdaac.usgs.gov/products/astgtmv003/), was used for elevation extraction with a spatial 
resolution of 30 m. 


2.3 Methods 


2.3.1 Snow cover parameters 

We analyzed the spatiotemporal variation characteristics of snow cover on the Mongolian Plateau 
with several snow cover parameters, such as the SCS, SCF, SCD, SCOD, and SCED. We defined 
the period from 1 September of the current year to 30 April of the next year as an SCS; for 
example, period from 1 September 2000 to 30 April 2001 was the SCS of 2001, which can be 
relative to the time range of grassland phenology parameters. MODIS snow products C6 does not 
provide SCF dataset; thus, we calculated the SCF of the study area according to the equation in 
the MODIS snow products C6 user guide (Riggs et al., 2015). The formula is as follows: 


SCF = (—0.01+ (1.45 x NDSI)) x10096, (1) 


where NDSI is Normalized Difference Snow Index. The total number of snow days in the study 
area in an SCS is defined as the SCD. The SCOD is the day of year (DOY) when snow first 
appeared in any pixel and the following five consecutive days were covered by snow during an 
SCS. The SCED is the DOY when any pixel has snow for the last time, and the pixel was covered 
by snow for five consecutive days before. 

The maximum value composite method was used to calculate the monthly SCF, and then the 
average value method was used to calculate the multi-year average of each snow cover parameter. 

We classified the SCOD and SCED according to the month scale, i.e., 1-31 d corresponds to 
January, 32—59 d corresponds to February, 60—90 d corresponds to March, 91—120 d corresponds 
to April, 244—272 d corresponds to September, 273—304 d corresponds to October, 305—335 d 
corresponds to November, and 336—365 d corresponds to December. It is worth noting that SCOD 
or SCED might be too late or too early, causing discontinuities in time. Therefore, SCOD later 
than 365 DOY is uniformly classified as next year, while SCED earlier than 1 DOY is uniformly 
classified as last year. 
2.3.2 Grassland phenology 


The phenology parameters were calculated by using the extreme value of the logistic curve 
curvature of the cumulative NDVI, which fitted the cumulative NDVI value in the year to simulate 
the NDVI change curve (Bao et al., 2017) using the logistic function. The formula is as follows: 


C 
t) 2 ———--* d, 2 
y(t) lie (2) 


where t is the Julian day; y(t) is the cumulative NDVI for logistic fitting corresponding to t; a and b 
are fitting parameters; c is the maximum accumulated NDVI; and d is the initial background NDVI 
value. Then, the logistic fitting curve curvature (K) of the accumulated NDVI was calculated 
(Zhang et al., 2003; Hou et al., 2014b; Bao et al., 2017). Finally, the SOS and the end of the 
growing season (EOS) of grassland phenology were extracted by pixel though the curvature 
extreme value (Bao et al., 2017). The formula is as follows: 
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where a is the value of the accumulated NDVI; s stands for time; z=e@*9; and a, b, c, and t are the 
same as in Equation 2. The SOS and EOS determine the DOY when the cumulative curve curvature 
reaches its maximum and minimum values. LOS is the time between the beginning and the end of 
the vegetation growth season. According to the actual situation of the plateau and combined with the 
results of previous study (Bao et al., 2020), we believe that the region of NDVI<0.08 is the area 
without vegetation growth. Finally, we obtained the phenology parameter data from 2001 to 2018 
by using the average value method. 

2.3.3 Grey relation analysis (GRA) 


To quantify the effect of snow cover on grassland phenology, we calculated the grey relation grade 
(GRG) and determined the order of correlation between grassland phenology and snow cover 
parameters using the GRA (Li et al., 2016; Chang et al., 2018; He et al., 2018). The GRA establishes 
the correlation degree between the corresponding sequences according to the similarity degree of 
the geometric shape of the comparison sequence and the reference sequence curve; the closer the 
curve is, the greater the correlation degree is. In this paper, the snow cover parameter was set as a 
comparison sequence, X;=x;(k) (where i=1, 2, 3, 4, and k=1, 2, ..., n), and the phenology parameter 
was set as a reference sequence, Xo. To eliminate the impact of the dimension, we standardized the 
processing of the original snow cover and phenology data as follows: 


X,(k) - Xi 
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where X;iand X» are the means of the comparison sequence and reference sequence, respectively; 
and o; and oo are the standard deviations of the comparison sequence and reference sequence, 
respectively. The formula for calculating GRG is as follows: 

Ed Amin-- gA max l (6) 

' Ao, (t) + «^ max 

where &,,(k) is the grey relation coefficient; Ad((t) is the absolute value of the difference between the 
comparison sequence and the reference sequence; and Amin and Amax are the minimum and 
maximum values of Ad((t), respectively. a is the resolution coefficient, with the recommended value 
range of 0.1-0.5; the smaller the resolution coefficient, the better the discrimination effect. 
Therefore, this study set the resolution coefficient a=0.1 (Li et al., 2016). Finally, the mean of the 
grey relation coefficient was calculated, and the GRG was obtained. The formula is as follows: 
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where Z(k) is the grey relation coefficient; and r; is the GRG of the comparison sequence and the 
reference sequence. A larger value of the GRG indicates that the effect of the comparison sequence 
on the reference sequence is more obvious. 


2.3.4 Trend analysis 


The parameter change rates of snow cover and phenology parameters on the Mongolian Plateau 
from 2001 to 2018 were calculated using the unitary linear regression method, and the spatial 
variation trend of each parameter was studied. The formula is as follows: 
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where Ope denotes the change trend of snow cover and phenology parameters; n is the time series 
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length; i is the year sequence number; and X; is the value of each parameter in year i. Asiope>O 
indicates a growth trend, and Osiope<O indicates a decreasing trend. The results were tested using the t 
test, and the significance level was set as a=0.05. 


3 Results 


3.1 Spatiotemporal characteristics of snow cover on the Mongolian Plateau 


3.1.1 Spatial characteristics of snow cover 


Using MOD10A1 snow cover data, we obtained the spatial distribution information of the snow 
cover parameters and calculated the trend of snow cover parameter changes on the Mongolian 
Plateau from 2001 to 2018 (Fig. 2). Figure 2a shows that the average SCF of the Mongolian 
Plateau was 29.9996, and the spatial distribution exhibited a distribution of more in the north and 
less in the south. The high-value areas were concentrated in the Hangay Uul Mountains in the 
west, the Sayan Mountains in the northwest of the Mongolian Plateau, the Khentii Mountains in 
the north, and the Da Hinggan Mountains in the east, which are plateau and mountainous terrain, 
creating favourable conditions for snow cover accumulation, and the SCF can reach 75.0096 in 
local areas. The low-value areas were concentrated in the central and southern deserts of the 
Mongolian Plateau, the Gobi region and parts of the southeast. In addition, we found that there 
were obvious differences in the SCF between the eastern and western sides of the eastern Da 
Hinggan Mountains, and the SCF in the eastern forest area value was significantly smaller than 
that in the western grassland vegetation area. It is possible that the eastern region was affected by 
the forest canopy, which led to the underestimation of MODIS snow cover production. Figure 2e 
shows that the SCF in most areas of the Mongolian Plateau tended to decrease. In 73.5196 of the 
study area, the SCF tended to decrease and was mainly distributed in the southern and northern 
Khentii Mountains, eastern part of the Hulun Buir Plateau, and the Da Hinggan Mountains. Of 
these regions, 4.5396 showed a significant decrease and were mainly distributed in the southern 
part of the Altay Mountains. Twenty-six percent of the study area showed an increasing trend of 
the SCF, which were mainly distributed in the northwest of the Altay Mountains, the Hangay Uul 
Mountains, and the southeastern regions. Overall, the SCF of the Mongolian Plateau showed a 
downward trend with an annual change rate of —0.33% (+0.56%). 

Figure 2b indicates that the average SCD on the Mongolian Plateau was 57.96 d, and the spatial 
distribution was similar to that of the SCF; both of them showed a long-term SCS in the north and 
a short SCD in the south. The northern parts were stable snow cover regions, with the SCD above 
60 d, even reaching 120 d in parts of the northwest and northeast of the Mongolian Plateau. The 
central, southern and southeastern regions were unstable snow cover regions, with less than 60 d 
of the SCD. Figure 2f indicates that the SCD in most regions of the plateau was decreasing. A 
total of 70.73% of the regions showed a shortening trend, mainly concentrated in the north and 
south; 29.27% exhibited a lengthening trend, mainly distributed in the northwest and southeast. 
Overall, the SCD on the Mongolian Plateau showed a shortening trend with an annual change rate 
of —0.57 (+6.21) d. 

The SCOD over the plateau was very long and distributed between 305 and 365 DOY, showing 
early distribution characteristics in the north and late distribution in the southern regions (Fig. 2c). 
The areas with the early SCOD were distributed in the mountainous areas of the northwest, of 
which a small part of the SCOD were particularly early, between 244 and 272 DOY. The area had 
a high altitude, low temperature, and early time of occurrence. The SCOD in the northern and 
eastern regions was mainly focused approximately 305-335 DOY, while the SCOD in the central 
and southwestern regions was mainly concentrated in 336-365 DOY. In the deep inland of the 
above areas, changes in precipitation patterns had a substantial effect on snow cover. According to 
Figure 2g, the SCOD in 67.0096 of the Mongolian Plateau showed an advanced trend, mainly in 
the northwest, northeast and south, with an annual change rate of —5.33 (+5.17) DOY, and 33.00% 
of the area showed a delayed trend of the SCOD, which was irregular in the north and southeast, 
with an annual change rate of 2.59 (+2.85) DOY. 
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Fig. 2 Spatial distributions of the snow cover fraction (SCF; a), snow cover duration (SCD; b), snow cover 
onset date (SCOD; c), and snow cover end date (SCED; d), as well as the variations of the SCF (e), SCD (f), 
SCOD (g), and SCED (h) in the snow cover season (SCS) on the Mongolian Plateau from 2001 to 2018. DOY, 
day of year. The SCOD later than 365 DOY is uniformly classified as next year, while the SCED earlier than 1 
DOY is uniformly classified as last year. 
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The SCED on the Mongolian Plateau was equally long and distributed between 16 and 120 
DOY, showing the late distribution characteristic of the south and early distribution characteristic 
in the north (Fig. 2d). The early SCED in the south and southeast are focused approximately 1-31 
DOY. Compared with other areas, these regions had high temperatures, which affected snow 
cover. The SCED in the northern region was late and focused approximately 60-120 DOY; this 
region was in the high latitudes with low temperatures and strong snow protection. Figure 2h 
shows that 72.8496 of the regions on the Mongolian Plateau showed an advancing trend of the 
SCED, mainly in the south and east, with an annual change rate of —5.74 (+5.58) DOY. A total of 
27.1696 of the area showed a delayed trend, mainly in the northwest, with an annual change rate 
of 2.83 (+3.20) DOY. 

3.1.2 Interannual characteristics of snow cover 

The variation in snow parameters during different SCSs was not stable and showed a very intense 
wave. As shown in Figure 3, the SCFs fluctuated between 27.0196 and 53.9696 over the last 18 a; 
the minimum SCF was 27.1196 in 2018, and the maximum was 53.8496 in 2002 (Fig. 3a). The 
SCD fluctuated between 33 and 80 d from 2001 to 2018; the minimum SCD was 33 d in 2018, 
and the maximum was 80 d in 2012 (Fig. 3b). The SCOD fluctuated between 317 and 343 DOY 
during the period of 2001—2018; the earliest SCOD was 318 DOY in 2012, and the latest was 343 
DOY in 2005 (Fig. 3c). The SCED fluctuated between 49 and 66 DOY from 2001 to 2018; the 
earliest ending day was 49 DOY in 2005, and the latest was 66 DOY in 2012 (Fig. 3d). Overall, 
the SCF and SCD showed a nonsignificant downward trend, and the SCOD and SCED showed a 
nonsignificant postponement trend. 
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Fig.3 Temporal changes in the SCF (a), SCD (b), SCOD (c), and SCD (d) in the SCS on the Mongolian Plateau 
from 2001 to 2018 


3.2 Spatiotemporal characteristics of vegetation phenology on the Mongolian Plateau 


3.2.1 Spatial characteristics of vegetation phenology 

Based on MOD13A1 NDVI data, we obtained the spatial distribution and change trend in 
grassland phenology on the Mongolian Plateau from 2001 to 2018 by using the extreme value of 
logistic curve curvature of the cumulative NDVI (Fig. 4). The spatial distribution of grassland 
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phenology was similar to that of snow cover, showing obvious differences between the north and 
the south. Figure 4a shows the spatial distribution of the average SOS of the Mongolian Plateau 
from 2001 to 2018. The SOS was concentrated in the range of 100-140 DOY, showing the spatial 
distribution characteristics of a delay from the south to the north, with a strong spatial 
heterogeneity. The SOS occurred early in the central and southern Gobi Desert areas, mostly on 
91-110 DOY. In contrast, the SOS in the northern region occurred late, mostly on 121-140 DOY, 
including the western Hangay Uul Mountains, the Sayan Mountains, and the eastern Da Hinggan 
Mountains; in some parts of the Hangay Uul Mountains and the Da Hinggan Mountains, the SOS 
occurred after 140 DOY. According to Figure 4d, more than half of the regions (51.9096) on the 
Mongolian Plateau showed different degrees of advancement of the SOS, which were mainly 
distributed in the northwest of the Altay Mountains and the Sayan Mountains, the northern 
Hangay Uul Mountains, and the eastern Da Hinggan Mountains, with an annual change rate of 
—0.39 (+0.37) DOY. Forty-eight percent of the regions showed different degrees of delay in the 
SOS, mainly in the central region, with an annual change rate of 0.75 (+7.47) DOY. 

Figure 4b shows the average spatial distribution of the LOS on the Mongolian Plateau from 
2001 to 2018. The LOS of the Mongolian Plateau was concentrated during 121—180 d, with the 
spatial distribution characteristics of a short growth season in the north and a long growth season 
in the south. The Altay Mountains, the Hangay Uul Mountains, and the Sayan Mountains in the 
northwest had the LOS of 100—120 d, but the mountain areas had an LOS of 141—160 d, which 
was approximately 20 d longer than that in the above areas. The LOS in other areas in the north 
was between 101 and 180 d, which was mainly manifested by the short LOS in the Hulun Buir 
Plateau of between 101 and 120 d. The LOS in the Khentii Mountains and the eastern part of the 
Da Hinggan Mountains was long, ranging from 161 to 180 d. The LOS in the central and southern 
regions was relatively long, exceeding 161—180 d, and the LOS in some areas even exceeded 181 
d. Figure 4e shows that the LOS in most areas of the Mongolian Plateau showed an extended 
trend. A total of 59.7096 of the areas showed an increasing trend of the LOS, which was mainly 
distributed in the north and south. The southern side of the Altay Mountains in the northwest, the 
Khentii Mountains in the north, and the northern side of the Da Hinggan Mountains in the east 
showed a significant extended trend; 40.3096 of the areas showed a shortened trend of the LOS, 
mainly distributed in the central and eastern Hulun Buir Plateau. Overall, the LOS on the 
Mongolian Plateau showed a prolongation trend, with an annual change rate of 0.19 (+1.29) d. 

Figure 4c displays the average spatial distribution of the EOS on the Mongolian Plateau from 
2001 to 2018. The EOS on the Mongolian Plateau was concentrated from 261 to 280 DOY, 
showing the spatial distribution characteristics of early in the north and late in the south. The EOS 
appeared early in the northern region, mostly between 256 and 270 DOY, especially in the Altay 
Mountains, the Hangay Uul Mountains, and the Sayan Mountains. In a small area of the eastern 
Da Hinggan Mountains in the northwest, the EOS occurred very early, between 185 and 255 
DOY. In contrast, the EOS occurred late in the central and southern regions, between 276 and 285 
DOY. Figure 4f shows that the EOS in most areas of the Mongolian Plateau was delayed. Of all 
these regions, 64.3196 had delayed EOS, mainly distributed in the central Gobi Desert, the 
northern Khentii Mountains, and part of the eastern Da Hinggan Mountains. In 35.6996 of the 
regions, the EOS showed an advancing trend, mainly distributed in the northwestern Altay 
Mountains, the Sayan Mountains, the southern Ordos Plateau, the Ulanqab Plateau, and the 
eastern Hulun Buir Plateau. Overall, the EOS on the Mongolia Plateau showed a deferred trend, 
with an annual change rate of 0.31 (+5.17) d. 

In addition, similar to the spatial distribution of snow cover, there were significant differences 
in the SOS, LOS, and EOS between the eastern and western Da Hinggan Mountains. The SOS in 
the eastern forest-covered area was earlier than that in the western inland grassland area, the EOS 
was later than that in the western inland grassland area, and the LOS was longer than that in the 
western inland grassland area. The Da Hinggan Mountains may block the abundant water vapour 
from entering the inland Pacific Ocean, which makes the temperature and precipitation conditions 
in the inland grassland area unable to meet the requirements of grassland phenology and gives the 
phenology in the Da Hinggan Mountains an obvious spatial heterogeneity. 
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Fig. 4 Spatial distributions of the start of growing season (SOS; a), the length of growing season (LOS; b), and 
the end of the growing season (EOS; c), as well as the changes in SOS (d), LOS (e), and EOS (f) on the 
Mongolian Plateau from 2001 to 2018 


3.2.2 Interannual characteristics of vegetation phenology 


In the time series, the change of grassland phenology parameters showed a dramatic inter-annual 
fluctuation (Fig. 5). In the past 18 a, the SOS has ranged from 112 to 130 DOY, with no 
significant changes. The earliest SOS was 112 DOY in 2009, and the latest SOS was 130 DOY in 
2003. The EOS was between 263 and 279 DOY, showing a non-significant postponement trend, 
among which the earliest EOS was 263 DOY in 2009, and the latest EOS was 279 DOY in 2018. 
The LOS ranged from 140 to 158 d, showing a non-significant prolongation trend, with the 
shortest value being 140 d in 2003 and the longest value being 158 d in 2007. 
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Fig. 5 Temporal changes in the grassland phenology of SOS (a), EOS (b), and LOS (c) on the Mongolian 
Plateau from 2001 to 2018 


3.3 Impact of snow cover changes on grassland phenology 


Table 1 shows the correlation coefficient between the snow cover parameters of the Mongolian 
Plateau and the SOS of four grassland vegetation types from 2001 to 2018. The SCF, SCD, and 
SCED were significantly positively correlated with the SOS of the four types of grassland, among 
which the SCF had the closest correlation with meadow steppe, the SCD had the closest 
correlation with desert steppe, and the SCED had the closet correlation with typical steppe. The 
SCOD was significantly negatively correlated with the SOS of the four types of grassland and 
was most strongly correlated with the SOS in the alpine steppe. This result means that with the 
increase in the SCF, the extension of the SCD, and the postponement of the SCED, the SOS of the 
four types of grassland vegetation will be postponed, among which the increase in the SCF will 
delay the SOS of the meadow steppe the most, the extension of the SCD will delay the SOS of the 
desert steppe the most, and the delay in the SCED will delay the SOS of the typical steppe the 
most; however, the postponement of the SCOD will advance the SOS of the four types of 
grassland vegetation, among which the SOS of the alpine steppe will be the most advanced. 


Table 1 Pearson correlation coefficient between the snow cover parameters and the start of growing season 
(SOS) of the four types of grassland on the Mongolian Plateau from 2001 to 2018 


Grassland type SCF SCD SCOD SCED 
Meadow steppe 0.649" 0.489" —0.300"* 0.368" 
Typical steppe 0.586" 0.632" -0.336" 0.540" 
Desert steppe 0.585* 0.647* -0.267* 0.497* 
Alpine steppe 0.5677 0.585" —0.393™ 0.522™* 


Note: SCF, snow cover fraction; SCD, snow cover duration; SCOD, snow cover onset date; SCED, snow cover end date. **, significant 
correlation at P«0.01 level. 


The GRG value can indicate the influential intensity of different snow cover parameters on the 
SOS of grassland vegetation. Overall, the SCD had the greatest impact on the SOS of the four 
grassland vegetation types on the Mongolian Plateau (Fig. 6). For the meadow steppe, the mean 
GRG of the SCF was the largest (0.7152); for the typical steppe, desert steppe, and alpine steppe, 
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the mean GRG of the SCD reached its maximum value (0.7152, 0.7100, and 0.6891, 
respectively), which shows that the SCF had a significant effect on the SOS of the meadow 
steppe, and the SCD had the greatest impact on the SOS of the typical steppe, desert steppe, and 
alpine steppe. The SCOD and mean GRGs of the four grassland vegetation species reached their 
minimum values (0.5774, 0.5639, 0.5441, and 0.5121, respectively), which indicated that the 
SCOD had the least impact on their SOS. However, the residual snow parameters had different 
effects on the SOS of different grassland vegetation types. For typical steppe and alpine steppe, 
the SCED affected the vegetation SOS to a lesser extent than the SCD, and the SCF had the 
second greatest effect. The impact of the SCF on the SOS of the desert steppe was second only to 
that of the SCD, followed by the SCED. For the SOS of the meadow steppe, the impact of the 
SCF was second only to that of the SCD, followed by the SCED. 


SCF 


SCD "----  SCED 
SCF 


Meadow steppe Typical steppe Desert steppe Alpine steppe 


Fig. 6 Grey relation grade (GRG) values between the snow cover parameters and the SOS for the four grassland 
vegetation types on the Mongolian Plateau from 2001 to 2018 


Tables 2 and 3 present the correlations of the snow cover parameters with the EOS and LOS of 
the four grassland vegetation types from 2001 to 2018. The correlation of the three snow cover 
parameters (SCF, SCD, and SCED) with EOS and LOS is just opposite to that with SOS. The 
specific manifestation is that the SCF, SCD, and SCED were significantly negatively related to 
the EOS of the four grassland types, while the SCOD and the EOS showed a significant positive 
correlation. The Pearson correlation coefficient between the snow cover parameters and the EOS 
of the alpine steppe reached a maximum value. This finding indicates that snow was closely 
related to the EOS of the alpine steppe. The increase in the SCF, the extension of the SCD, and 
the postponement of the SCED will lead to the advancement of the EOS in the alpine steppe, 
while the postponement of the SCOD will postpone the EOS in the alpine steppe. In addition, 
although the snow cover parameters had significant correlations with the meadow steppe, typical 
steppe, and desert steppe, the correlation coefficients were small, indicating that the EOS changes 
in these three grassland vegetation types showed a weak correlation with snow cover. 


Table 2 Pearson correlation coefficient between the snow cover parameters and the end of the growing season 
(EOS) of the four grassland vegetation types on the Mongolian Plateau from 2001 to 2018 


Grassland type SCF SCD SCOD SCED 

Meadow steppe —0.298"* —0.138" 0.0817 -0.112* 
Typical steppe -0.275" -0.275" 0.202" -0.263" 
Desert steppe -0.193* -0.201* 0.150* -0.133* 
Alpine steppe -0.484* —0.489'* 0.425" -0.402* 


Note: **, significant correlation at P«0.01 level. 


JOURNAL OF ARID LAND 2021 Vol. 13 No. 4 


Table 3 Pearson correlation coefficient between the snow cover parameters and the length of the growing 
season (LOS) of the four grassland vegetation types on the Mongolian Plateau from 2001 to 2018 


Grassland type SCF SCD SCOD SCED 

Meadow steppe —0.626™ -0.418* 0.254* —0.319* 
Typical steppe —0.615™ —0.650** 0.381" 0.573" 
Desert steppe —0.624"* —0.683" 0.321” -0.511" 
Alpine steppe —0.635* —0.650** 0.487* —0.562* 


Note: **, significant correlation at P<0.01 level. 


The correlation between the snow cover parameters and the LOS was similar to that of the EOS 
(Table 3). The SCF, SCD, and SCED were significantly negatively related to the LOS of the four 
grassland vegetation types, while the SCOD had a significant positive correlation with the LOS, 
among which the SCF had the closest correlation with the LOS of the alpine steppe, the SCD had 
the closest correlation with the LOS of the desert steppe, and the SCED had the closest 
correlation with the LOS of the typical steppe. These findings indicate that with the increase in 
the SCF, the extension of the SCD and the extrapolation of the SCED will occur later, and the 
LOS of all four types of grassland will be shortened. The increase in the SCF will shorten the 
LOS of the alpine steppe the most, the extension of the SCD will shorten the LOS of the desert 
steppe the most, and the delay in the SCED will shorten the LOS of the typical steppe the most. 
The correlation between the SCOD and LOS was negative for all grassland vegetation types, 
which means that the delay in the SCOD will lead to the extension of the LOS, among which the 
LOS of the alpine steppe will be the most prolonged. 

Overall, the SCOD and SCED had a greater impact on the EOS and LOS of the four grassland 
vegetation types on the Mongolian Plateau, while the SCF and SCD exhibited a smaller impact. 
Figure 7 shows that the order of GRG values of the EOS, LOS, and snow cover parameters of 
different grassland vegetation types was very similar, indicating that the influences of different 
snow cover parameters on the EOS and LOS of different grassland vegetation types are roughly 
the same. For the four types of grassland vegetation, the SCOD had the greatest impact on the 
EOS and LOS, followed by the SCED. For the meadow steppe, the SCF had the weakest effect on 
the EOS and LOS. For the other three grassland vegetation types, the SCD had the weakest 
influence on the EOS and LOS, which may be due to the lag effects of both the season and 
vegetation. 
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Fig. 7 GRG values between the snow cover parameters and the EOS (a) and between the snow cover 
parameters and the LOS (b) of the four grassland vegetation types on the Mongolian Plateau from 2001 to 2018 
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4 Discussion 


The Mongolian Plateau is one of the most snow-rich regions in Central Asia because of its cold 
and long winters. This study used MODIS snow cover data to analyze the snow cover situation on 
the Mongolian Plateau during 2001-2018. The results shown in Figure 2a-d indicated that the 
SCF was high, the SCD was long, the SCOD was early, and the SCED was late on the northern 
plateau, and the snow cover situation in the south was the opposite of that in the north (Sa et al., 
2012; Han et al., 2014). This difference may exist because the north is affected by Mongolian 
high pressure during the snow season and temperatures in the north are considerably lower than 
those in other regions; additionally, water vapor from the Arctic Ocean, the Black Sea, the Yellow 
Sea, and the Bohai Sea can easily reach this area to form precipitation. Further, the area is located 
at high latitudes, the solar radiation energy received during the snow season is relatively limited, 
and the snow cover reflects most of the solar radiation that travels to the surface, resulting in low 
surface temperatures. In the Gobi and desert areas of the central and southern plateau, it is 
difficult for water vapor to reach this area to form precipitation because of its inland location; the 
temperature in the above regions is higher than that in other areas in the same time range, which 
is not conducive to the formation of snow cover. In addition, both the sublimation and 
evaporation of snow cover can cause the surface snow cover to decrease (Zhang et al., 2016). 
These factors result in the small SCF, the short SCD, the late SCOD, and the early SCED in the 
central and southern plateau, and thus, there is a significant difference in the snow distribution 
from that in the north. 

In this study, the grassland phenology of the plateau was analyzed between 2001 and 2018 
using MODIS NDVI data, and the results show that the spatial heterogeneity of the grassland 
phenology on the plateau was also obvious due to the topography, altitude, latitude, climate 
difference, and human activities. At the spatial scale, the grassland phenology (late SOS, early 
EOS, and short LOS) in the northern regions (Fig. 4a—c) was basically consistent with the results 
of Bao et al. (2017). This distribution pattern may be due to the high latitude and elevation in the 
north, which will result in lower temperatures, finally affecting vegetation growth. Although the 
southern regions were dry with less rainfall and high evapotranspiration, they showed the 
distribution of early SOS, late EOS, and long LOS compared with the northern regions, which 
may be due to the impact of human activities. Animal husbandry in this area (Tong et al., 2018) is 
the main economicactivity, and an earlier SOS can provide the large amount of forage needed for 
spring animal husbandry as early as possible. The extension of the grassland growth period is of 
great importance to herdsmen and local economies (Miao et al., 2017). Considering the natural 
factors, the rising temperature and precipitation or spring snowmelt water due to global warming 
are the two main factors leading to an early SOS, a late EOS, and a long LOS of vegetation in the 
south (Gong et al., 2015). During the SCS, snow can preserve the surface heat and protect crops 
(mainly winter wheat) from cold air (Chen et al., 2015). An earlier SCED will advance the 
snowmelt date, which also will lead to significant changes in spring snowmelt time and runoff 
and increase the possibility of catastrophic events, such as spring flood and summer drought. 

Affected by global warming, the snow cover area in the Northern Hemisphere is undergoing a 
decreasing trend (IPCC, 2008; Brown and Mote, 2009; Brown and Robinson, 2011). In recent 
years, many studies have shown that the warming temperatures have reduced SCDs and snow 
depth (Gao et al., 2012; Guo and Li, 2015; Ke et al., 2016; Xu et al., 2017), which is sure to affect 
the soil moisture in grassland areas and spring vegetation green up, as well as interfering with 
vegetation growth and affecting terrestrial ecosystems (Han et al., 2015). Therefore, we analyzed 
the snow cover and grassland phenology on the Mongolian Plateau and found that the SCF and 
SCD were significantly positively correlated with the SOS, and the correlation coefficient was 
large (Table 1). This conclusion is consistent with those of previous studies (Xie et al., 2017; 
Wang et al., 2018a). Figure 6 shows that the GRG values of the SCF and SCD to the SOS of the 
grassland vegetation are greater than those of the SCOD and SCED. These results indicate that 
the SCF and SCD are the main influencing factors affecting grassland vegetation green up on the 
plateau. This situation exists because in the cold winter, the snow that covers the surface will 


JOURNAL OF ARID LAND 2021 Vol. 13 No. 4 


insulate the surface vegetation from the air, protect the vegetation from dying due to strong winds 
and cold temperatures, and play a warming role for the plants (Chen et al., 2015; Qiao and Wang, 
2019). In addition, snowmelt water becomes an important source of water for plants when the 
snow melts in spring (Guo and Li, 2015; Wang et al., 2018a). In this study, the SCF was used as a 
parameter to characterize the snow cover range on the surface, and the SCD was used as a 
parameter to illustrate the SCD on the surface. These values can well indicate the growth 
conditions of grassland vegetation in spring under a natural state. 

Similarly, we determined the relationship between the snow cover parameters and the EOS and 
LOS of the four grassland vegetation types and found that although the snow cover parameters 
were significantly related to the EOS and LOS, the correlation coefficients were small (Tables 2 
and 3), which indicates that although snow cover had a certain influence on the EOS and LOS of 
grassland vegetation, it was not the main factor (Sha et al., 2016). The order of GRG values in the 
snow cover parameters and the EOS as well as the LOS shows that the SCOD is the main snow 
cover factor affecting the EOS and LOS of grassland vegetation. This result may arise because the 
early and late SCOD may affect the physiological indices of grassland vegetation (leaf area size, 
root development degree, and so on), thus affecting the EOS and LOS (Wang et al., 2018a). For 
example, an earlier SCOD will result in frostbite on the surface parts of grassland vegetation, 
causing irreversible damage to vegetation and affecting turn green and subsequent growth and 
development of plants in the second year. At the same time, the SCOD can indirectly reflect the 
temperature and precipitation on the plateau in autumn and winter, and the temperature and 
precipitation are important factors affecting the growth of grassland vegetation. The SCED is the 
second snow cover factor affecting the EOS and LOS of grassland vegetation, which may affect 
the EOS and LOS by influencing the SOS of grassland vegetation (Wang et al., 2018a). In 
addition, the change in winter snow cover can change the soil water and heat conditions, thus 
affecting the nutrients needed for vegetation growth and the length of the vegetation growth 
period (Wang et al, 2018a). In future study, the processes of snow cover, snow melt, and 
vegetation growth will be fully monitored, and the mechanism of snow cover that affects 
vegetation growth on the Mongolian Plateau grasslands will be clarified. 


5 Conclusions 


In this study, the spatiotemporal distributions and change trends of snow cover parameters and 
grassland phenology on the Mongolian Plateau were analyzed using MOD10A1 snow cover data 
and MOD13A1 NDVI data from 2000 to 2018. We used a correlation analysis and GRA to 
investigate the effects of snow cover parameters on grassland phenology on the Mongolian 
Plateau. The conclusions are outlined as follows: 

(1) Winter snow cover exhibited a significant different spatial distribution on the Mongolian 
Plateau from 2001 to 2018. Compared with the southern part of the Mongolian Plateau, the 
northern part had a higher SCF, a longer SCD, an earlier SCOD, and a later SCED. The SCF and 
SCD showed non-significant downward trends at an annual rate of —0.33% (30.5696) and —0.57 
(+6.21) d, respectively. In total, the SCOD and SCED exhibited non-significant postponement 
trends in more than half of the Mongolian Plateau. The annual change rates were —2.63 (+9.49) 
and —4.01 (£11.27) DOY, respectively. 

(2) From 2001 to 2018, the average SOS and EOS occurred on 121 DOY and 271 DOY, 
respectively. They showed an opposite spatial distribution, i.e., the SOS was delayed from the 
south to the north, while the EOS delayed from the north to the south. The average LOS lasted for 
150 d; it showed a spatial distribution of delayed from the north to the south. The SCF, SCD, and 
SCED showed significant positive correlations with the SOS while significant negative 
correlations with the EOS and LOS for the meadow steppe, typical steppe, desert steppe, and 
alpine steppe. By contrast, the SCOD was significantly negatively correlated with the SOS and 
significantly positively correlated with the EOS and LOS. 

(3) The result of GRA shows that the SCF had the greatest impact on the SOS of the meadow 
steppe, the SCD had the greatest impact on SOS of the typical steppe, desert steppe, and alpine 
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steppe. The SCOD had the greatest influence on the EOS and LOS of the four grassland types. 

This study analyzed the changes in snow cover and its impact on vegetation phenology on the 
Mongolian Plateau in recent 18 years, which can provide a basic support for coping with future 
climate change and vegetation growth in arid and semi-arid grassland ecosystem. In the next step, 
we will use land surface model to in-depth study the mechanism of snow cover change on 
vegetation phenology combined with multi-source remote sensing data, such as snow depth, snow 
water equivalent, soil temperature, soil water content, etc. 
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